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Quantum Gravity effect on the Quark-Gluon Plasma 
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The Generalized Uncertainty Principle (GUP), which has been predicted by various theories of 
quantum gravity near the Planck scale is implemented on deriving the thermodynamics of ideal 
Quark-Gluon Plasma (QGP) consisting of two massless quark flavors at the hadron-QGP phase 
equilibrium and at a vanishing chemical potential. The effective degrees of freedom and MIT bag 
pressure are utilized to distinguish between the hadronic and partonic phases. We find that GUP 
makes a non- negligible contribution to all thermodynamic quantities, especially at high tempera- 
tures. The asymptotic behavior of corresponding QGP thermodynamic quantities characterized by 
the Stephan-Boltzmann limit would be approached, when the GUP approach is taken into consid- 
eration. 
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I. INTRODUCTION 

Various theories of quantum gravity predict essential modifications in the Heiscnbcrg's uncertainty principle near 
the Planck scale. In this paper, we utilize the proposed generalized uncertainty principle (GUP), which proved 
compatible with string theory, doubly special relativity and black hole physics. It seems that this approach 
accordingly modifies almost all mechanical Hamiltonians. Therefore, it can be implemented on studying the 
thermodynamics. 

The existence of a minimal length is one of the most interesting predictions of some approaches related to 
quantum gravity such as string theory as well as black hole physics. This is a consequence of string theory 
since strings can not interact at distances smaller than their size which yields GUP Q. Black hole (BH) physics 
suggests that the uncertainty relation should be modified near the Planck energy scale because of measuring 
the photons emitted from the black hole suffers from two major errors. The first one is the error by Heisenbcrg 
classical analysis and the second one is because the black hole mass varies during the emission process and the 
radius of the horizon changes accordingly [H-Q • 

Recently, a new model of GUP was proposed It predicts a maximum observable momentum and a 

minimal measurable length. Accordingly, = [Pi 7 Pj] = (via the Jacobi identity) results in. 



[Xi,Pj\ 



ill 



PiPj 
P 



+ a 2 {p 2 8 l] + 3piPj) 



(1) 



where a = q.q/M p c — a^l v fh and M p c 2 stand for Planck energy. M p and l p is Planck mass and length, 
respectively ao sets on the upper and lower bounds to a. Apparently, Eqs. ([I]) imply the existence of a 
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minimum measurable length and a maximum measurable momentum 

w a Q £ P , (2) 
a ~ M p c fo\ 

&Pmax ~ ., (3) 

where Ax > Ax rn i„ and Ap < Ap max . Accordingly, for a particle having a distant origin and an energy scale 
comparable to the Planck's one, the momentum would be a subject of a modification [6|-|8| 

Pi = Poi (1 - ap + 2a 2 pl) , (4) 

where Xi = xoi and poj satisfy the canonical commutation relations [xoi,poj] = $ij and simultaneously fulfil 
Eq. (TT|). Here, pot can be interpreted as the momentum at low energies (having the standard representation in 
position space, i.e. poi — — ihd/dxoi) and pi as that at high energies. 

The proposed GUP is assuming that the space is discrete, and that all measurable lengths are quantized 
in units of a fundamental minimum and measurable length. The latter can be as short as the Planck length 
0| ■ I n order to support the idea of this procedure, we can mention that similar quantization of the length 
(spatial dimensions) has been studied in context of loop quantum gravity Q. Furthermore, it has been suggested 
recently [Tl| that the GUP implications can be measured directly in quantum optics lab which seems to confirm 
the theoretical predictions EH 

Since the GUP apparently modifies the fundamental commutator bracket between position and momentum 
operators, then it is natural to expect that this would result in considerable modifications in the Hamiltonian. 
Furthermore, it would affect a host of quantum phenomena, as well. It is important to make a quantitative 
study of these effects. In a series of earlier papers, the effects of GUP was investigated on atomic and condensed 
matter systems @, H34I3I . on the weak equivalence principle (WEP), and on the Liouville theorem (LT) in 
statistical mechanics |15| . For instance, it has been found that GUP can potentially explain the small observed 
violations of the WEP in neutron interferometry experiments [l6j. Also, it can predict the existence of a 
modified invariant phase space which is relevant to the Liouville theorem. 

In this paper, we present a study for the impact of the GUP on QGP consisting of two massless quark flavors 
at hadron-QGP phase equilibrium and at a vanishing chemical potential. Wc calculate the corrections to various 
thermodynamic quantities, like pressure. The effective degrees of freedom and MIT bag pressure are utilized in 
order to distinguish between the hadronic and partonic phases. 

This paper is organized as follows. In section [Til the thermodynamics of QGP consisting of two massless 
quark flavors is briefly reviewed. In section IIII1 the effects GUP [H, Q on QGP thermodynamics is outlined. 
The results and discussions are given section IPV1 



II. THERMODYNAMICS OF QUARK-GLUON PLASMA 

In this section, we briefly review the thermodynamics of QGP [13]. At finite temperature T and chemical 
potential fj,, the grand-canonical partition function zb for non-interacting massive bosons with g internal degrees 
of freedom is given as 



In z B = \ 



^exp 



' T 



1 'J 



J[l-exp (E(k)-fi)Y 



(5) 



where k is the momentum of the particle and / is the occupation number for each quantum state with energy 
E(k) = \J k' 2 + to 2 with mass m. Here the infinite product is taken for all possible momentum states. For 
simplicity we consider a massless pion gas in a vanishing chemical potential. Then the grand potential reads 



CI 

v= 9 



d 3 k 



Tin 



1 — exp 



-E(k) 



T 



= -g—T 



(6) 



Therefore, pressure P, energy density e, and entropy density S of hadronic state can be deduced 

Ph = 9^T\ (7) 

e H = 3 5 .^T 4 , (8) 

S H = 4g.^T 3 . (9) 
Using MIT bag model [HI, then the QGP thermodynamic quantities accordingly read 

Pqgp = g QG p^T 4 -B, (10) 

sqgp = 3g QGP ^T 4 + B, (11) 

Sqgp = 4 3QGP ^T 3 , (12) 

where B is the bag pressure. Based on Gibbs condition [l9[ , the critical point is obtained at the phase equilib- 
rium, which is satisfied when Ph(T c ) = Pqgp(T c ). Then, the bag pressure can given as 

B = (g QG p - 9^ T c (13) 



III. GUP AND QGP 



For a particle of mass M having a distant origin and an energy comparable to the Planck scale, the momentum 
would be a subject of a tiny modification and so that the dispersion relation would too @, H, [To|, 0] 



E{k) z = k A c 2 (l - 2ak) + MV, 



(14) 



where c and M p i are the speed of light as introduced by Lorentz and implemented in special relativity and 
Planck's mass, respectively. For simplicity we use natural units in which h = c = 1 and will consider a masslcss 
pion gas 

E(k) = k(l - 2ak) 1/2 (15) 
For large volume, the sum over all states of single particle can be rewritten in terms of an integral [TBj 
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(16) 



Therefore, the partition function, Eq. ([B]), reads 

In 
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(18) 



It is obvious that the first term in Eq. (fl~8)) vanishes. Thus 



In z B 



Vg f 00 fc 3 
2^ J 3(1 - akf 
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(19) 



Let x = —(1 — 2afc) 1 / 2 so that da; = — —rn^dk and 

T y ' T(l-2ak) 1 / 2 



1 

In Z B = 7: — o 



fc 3 



2?r 2 7 3(1 - afc) 3 e( K ) - 1 



(20) 



Apparently, as we are interested in the terms containing the first order of a, so x can be approximated as 
follows. 



x = £(l-2afc)V» 



fc(l - ak) 



and therefore 

k = xT + ak 2 w .tT + a(x 2 T 2 + 2axTk 2 + a 2 fc 4 ). 
When ignoring higher orders of a, then 

fc ss a,-T(l + axT). 

The partition function becomes 

ln^^f^ 3 T 3(1 + ^ T)3 ^ 



2tt 2 y 3 (1 - axTf e x -l' 



(21) 

(22) 
(23) 

(24) 



It is apparent that the integral contains Maclaurin series, which is a Taylor series expansion of a function about 




(1 + axTf 



1 + 6axT +[■■■ }a 2 x 2 T' 2 + [• • ■ }a 3 x 3 T 3 
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where [• • • ] stands for non-identical factors. When ignoring all terms cantoning a with order > 2, then 



In z B 



6tt 2 



00 , x 3 dx 
e x - 1 



6aT 4 



x dx 
e x - 1 



The partition function is related to the grand canonical potential, lnzs = — O/T, 

ft 

where/(0)^ is Bose and Fermi integrals 



-^T 4 [r(5)/(0) 5 -] - a -^T 5 [r(6)/(0) 6 -], 



^(y)n+i 



i 



i 



-dx, 



(25) 



(26) 



(27) 



(28) 



(29) 



r(n + 1) 7 (^ 2 + 2/ 2 ) 1/2 exp [(x 2 + y 2 ) 1 / 2 ] ± 1 
and r(n) = (n — 1)! is the gamma function. At y = 0, then 

/(0)± +1 = ^C(n)a± 

This formula is valid for integer (n > 2) with a+ = 1 — 2^ 1_n \a~ = 1. The Ricmann zeta function, £(n) = 

2 2 

I2?=i j { ~ n) with C(2) = ^r,C(3) = 1.202C(4) = ^r,C(5) = 1.037, etc. Substituting from Eq. {2SJ and Eq. {Ml) 

J 6 90 

into Eq. ([27]) we have, 
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-<?^T 4 - 24aAT 5 C(5) = -g^T 4 - ga.T 5 , 



90 



90 



(30) 
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where a\ = 24aC(5)/7r 2 . 

The pressure is directly related to the grand canonical potential, P = —fl/V. Then, in the hadronic phase 



P H = 9-T* + ga.T* 



(31) 



In a massless ideal hadronic gas, the energy density is related with the pressure by the equation of state, 
eg = 3 Pji , and the entropy density is given by the derivative of pressure with respect to temperature T 



s H = 3. 9 — T 4 + SgaiT 5 , 
S H = 4g^T 3 + 5g ai T 4 . 



(32) 
(33) 



Taking into consideration the relevant degrees of freedom, it is obvious that this set of equations, Eq. (|31|) . 
(|32|) and ([33]) . is valid in the hadronic state. In QGP state, the bag pressure should be inserted as done in Eq. 
(fTU)) and (fTTj) . From Gibbs phase equilibrium condition Ph{T c ) = Pqgp{T c ) [HI, the bag pressure B can be 
deduced 



— T c 4 + a x T c 5 = 



B 



90 - " gQGP - 9-k 

Then, at this value for the bag pressure, the thermodynamic pressure in QGP state reads 
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from which the normalized QGP energy density can be written as 



£qgp 



9 QGP 



90 



9QGP — 9n ( T t 



9QGP 



T 



gQGpaiT 



gQGP - 9tt ( T c 
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As given in Eq. ([12")) , the QGP entropy is not depending on the bag pressure, 

-j$T = A 9qgp^ + 5g QGP aiT 



(34) 



(35) 



(36) 



(37) 



IV. RESULTS AND DISCUSSIONS 



In Fig. [T] the QGP pressure is given as function of T. The degrees of freedom in the QGP state assuming 
Hf = 2 and n c = 3 are 9qgp = 37. In the pionic (hadronic) state, they are g^ — 3. The horizontal solid lines 
represent the Stephan-Boltzmann (SB) limits of the pressure. The solid curves give the results of the first terms, 
Eq. (|35|) . in which the GUP effects is not included. The second term, the term including a\, is drawn by dashed 
and dotted curves. These two curves are distinguished by the values of upper and lower bounds of a. As given 
in [l4| and Eq. ([T]), the first bound for ao is about ~ 10 17 , which would approximately gives a ~ 10~ 2 GeV -1 . 
The other bound of ao which is ~ 10 10 . This lower bound means that a ~ 1CP 9 GeV -1 . As discussed in [Toj . 
the exact bound on a can be obtained by comparing with observations and experiments (2p| . It seems that the 
gamma rays burst would allow us to set an upper value for the GUP-charactering parameter a (2lj . 

In right panel of Fig. [2l the lattice QCD results for p/T A for dynamic three (nj = 3) quarks are given in 
dependence on T and for different temporal dimensions of the lattice [22| . The lattice temporal dimension 
N T is directly related to the inverse temperature, 1/T = N T a, where a is the lattice spacing. We notice that 
for lattice QCD highest temperature T = 1000 MeV, the pressure, for instance, is almost 20% below the SB 
limit. Other lattice QCD simulations for different quark flavors seem to confirm that the thermodynamics of 
massive quarks is below SB 0, [Hj], left panel of Fig. [2] Various studies for the lattice QCD thermodynamics 
with different quark masses and flavors have been performed over the last 2 — 3 decades. For a recent review 
on the lattice QCD equation of state, the readers are referred to [25|, [2(|. In all these studies, the resulting 
thermodynamic quantities above T c arc found to be smaller than the SB limit of the quantity of interest. For a 
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lower GUP 



0.3 0.6 0.9 

T [GeV] 

Fig. 1: The pressure normalized to T 4 is given in dependence on T. The horizontal solid line represents the Stephan- 
Boltzmann (SB) limit. The solid curve gives the result of the first term, Eq. (|35p . in which the GUP effects is not included. 
The second term, the term including ai, is drawn by dashed and dotted curves. These two curves are distinguished by 
the upper and lower bounds of a. 



long time, this observation stands without a clear interpretation or at least the worldwide community was not 
unified on one single interpretation. At that time, the idea that hadrons above T c is conjectured to be dissolved 
into non-correlated almost independent constituents of quarks and gluons was widely distributed. The latter 
are believed to for an ideal QGP. With the RHIC discovery [27}, an equilibrated, but strongly coupled QGP 
is formed in the heavy-ion collisions. Therefore, the distinguish of its thermodynamics from SB turns to be 
eligible. It is widely believed that the SB limit can be approached when the stron g ru nning coupling a s nearly 
vanishes. For a recent experimental evaluation of a s , the readers are referred to [281 ] . In light of lattice QCD 
thermodynamics and a s (Q 2 ) and in case of neglecting the effects of quantum gravity, the partonic constituents 
would form an ideal gas at very high temperatures or energies. If the GUP approach is taken into consideration, 
the asymptotic behavior that is characterized by the SB limit would be reached than the earlier case. Based 
on this study, one would feel encouraged to implement GUP approach on QGP with dynamic thee quarks and 
having strong correlations and dissipativc properties [29r]37j . 




Fig. 2: Left panel shows the results the normalized pressure with the p4 action for different quark flavors [24|]. The data 
are obtained in a very coarse lattice (small temporal dimension N T — 4) with bare masses m q /T — 0.4 and m s /T — 1. 
Right panel gives the results for N T = 6,8 and 10 lattices for three quark flavors [S^]. Therefore, the SB limit reads 
PSb/T 4 « 5.21. It is indicated by an arrow. 
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Appendix A: Integration by parts for Eq. (|17p 



In order to solve Eq. ([TT|). let us assume that 



--(l-2 Qfe )V= 

u = in 1 — e J 



and 



dv = 



(1 - Oik) 



(Al) 



Then 



du = 



1 --(l-2afc) 1/2 

_T 

k 

- — (l-2afc)V2 

1 — e J 



1 — 3a/c 



(1 - 2afc) 1 /2 



dk, 



and 



3(1 -akf 



(A2) 



The partition function, Eq. (|17[) can be re-written as 



\nz B = --^5 V ut, lo° + 2~2 V 



(A3) 
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